Quantum logic operations and creation of entanglement in a 
scalable superconducting quantum computer with long-range 
constant interaction between qubits 

G. P. Berman^, A. R. Bishop^, D. I. Kamenev^, and A. Trombettoni^'^ 

Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 8754-5 and 
^ LN.F.M. and Dipartimento di Fisica, Universita di Parma, 
parco Area delle Scienze 7A Parma, L43100, Italy 

We consider a one-dimensional chain of many superconducting quantum interfer- 
ence devices (SQUIDs), serving as charge qubits. Each SQUID is coupled to its 
nearest neighbors through constant capacitances. We study the quantum logic op- 
erations and implementation of entanglement in this system. Arrays with two and 
three qubits are considered in detail. We show that the creation of entanglement 
with an arbitrary number of qubits can be implemented, without systematic errors, 
even when the coupling between qubits is not small. A relatively large coupling 
constant allows one to increase the clock speed of the quantum computer. We ana- 
lytically and numerically demonstrate the creation of the entanglement for this case, 
which can be a good test for the experimental implementation of a relatively sim- 
ple quantum protocol with many qubits. We discuss a possible application of our 
approach for implementing universal quantum logic for more complex algorithms by 
decreasing the coupling constant and, correspondingly, decreasing the clock speed. 
The errors introduced by the long-range interaction for the universal logic gates are 
estimated analytically and calculated numerically. Our results can be useful for ex- 
perimental implementation of quantum algorithms using controlled magnetic fluxes 
and gate voltages applied to the SQUIDs. The algorithms discussed in this paper 
can be implemented using already existing technologies in superconducting systems 
with constant inter-qubit coupling. 
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I. INTRODUCTION 

In the last few years, a large number of studies has been devoted to the realization of 
qubits using Josephson devices A single superconducting qubit can be realized using 
either charge or flux degrees of freedom. The flux Josephson two-state system is based on 
two quantum states carrying opposite persistent currents. Coherent time evolution between 
these states has been recently observed Q]- At present, however, no quantum oscillations 
between two coupled flux qubits has been reported. In charge Josephson qubits the rele- 
vant degree of freedom is the charge on superconducting grains. The coherent control of 
macroscopic quantum states in a single Cooper-pair box has been demonstrated j^, while 
the first observation of coherent quantum oscillations in two coupled charge qubits has been 
reportedj4l. Moreover, coherent oscillations have been observed in other superconducting 
devices Jy,Q]- While longer coherence times are desirable, these experiments show that the 
superconducting circuits are strong candidates for solid-state qubits. The next major step 
toward building a Josephson-j unction based quantum computer is to experimentally realize 
simple quantum algorithms, such as the creation of an entangled state involving more than 
two coupled qubits. 

A typical design of a Cooper pair box consists of a small superconducting island with 



n Cooper pair charges connected by a tunnel junction wit 



1 a Josephson coupling energy 



. A control gate voltage V is 



Ej and the capacitance Cj to a superconducting electrode 
coupled to the system via a gate capacitor Cg. 

The energy of the n Cooper pairs in the box is E = AEcn^ — 2enV , where Eq = e'^/2{Cg + 
Cj), and e is the electron charge. Considering later n as one of the canonical variables, and 
neglecting the constant term, we can represent E in the form E = 4Ec{n — , where 
Hg = CgV/2e is the gate charge (in units of 2e). When Eq ^ Ej, by choosing Ug close to 
the degeneracy point, ng = {2n + l)/2, only the states with n and n + 1 Cooper pairs play 
a role. In this case, the effective Hamiltonian of the two-state system can be written in the 
spin- 1 notation as 

Hi = -5"S" - 5^S^ (1) 

/1\ 

where the state with n Cooper pairs corresponds to the spin state and the state with 

\oJ 

n+1 Cooper pairs corresponds to the spin state ( 1 (see, for example, pj]); and S"^ are, 



3 



respectively, the z and x - components of spin-i operator; ~ 4£'c(l ~ 2^^) and ~ Ej 
are the effective magnetic fields which are controlled by the applied gate voltage and the 
magnetic flux. [See Eqs. (jHI) below.] 

The Hamiltonian of an array of coupled superconducting qubits can be written in the 
general form 

Af-l N-\ N-\ 

H = - E - E ^rsr + E f^s^sj, (2) 



i=0 i=0 i.J=0 



where a = x,y,z and is the number of qubits in the circuit. The explicit form of the 
coefficients f/^" in depends on the particular way in which the inter-qubit coupling is 
implemented. A variety of possible architectures to couple charge Josephson qubits has been 
proposed 0,0, [lO, 11, 12, [3, 13 • been suggested to realize a coupling of the type 



oc S^Sj using an inductor: all qubits were connected in parallel to an LC-oscillator which 
provided the two-qubit interaction. A possible limitation of this architecture is that the inter- 
qubit coupling term is valid only under the conditions that the phase conjugate to the total 
charge on the qubit capacitors fluctuates weakly, and that the eigenfrequency uq of the LC 
circuit is much larger than the typical frequencies of the qubit dynamics. Since scales with 
the number of qubits as 1/ ^/N, this limits the maximum number of qubits in the quantum 
register. A different proposed type of the inter-qubit coupling uses the Coulomb interaction 
between charges on the islands of the charge qubits [9J: this gives a coupling oc SfS^. A 
drawback of this approach is that one cannot switch the inter-qubit coupling in this scheme 
without introducing unwanted dephasing effects. Another proposed architecture employs 
two additional SQUIDs to connect each qubit: all Josephson charge qubits are coupled 
through a common superconducting inductance, and the resulting coupling is oc SfSJ 11 1. 
The same coupling of charge qubits can also be realized using mutual inductance 10 1. 
Finally, we mention that another experimental realization of a single superconducting charge 
qubit has been obtained using the current-biased Josephson junctions 0, Q]- A quantum 
computer architecture based on two capacitively coupled current-biased Josephson junctions 

3,Q- 



was discussed before [12, 



In this paper, we consider capacitively coupled SQUIDs, which currently are the only 
experimental setup that has provided an experimental detection of quantum oscillations in 
two coupled charge qubits In this architecture, the SQUIDs are connected via a constant 
capacitor, and the coupling is of the form oc S^S|. Varying the magnetic flux $ through the 
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FIG. 1: A schematic illustration of array of capacitively coupled SQUIDs. 



SQUIDs allows one to control the Josephson energy Ej, which is equivalent to controlling 
the effective magnetic field -B^ [Ij. Therefore, we will assume that it is possible to have 
the Josephson energies [i.e., coefficients in Eq.©] independently and locally variable in 
time. 

We note that experimentally the fast independent manipulation of gate voltages is pos- 
sible jisl . 16^ . At the same time, the fast independent control of magnetic fluxes through 
the different SQUIDs with the distance ~ 1 /im between them and the switching time 
~ 1 ns represents a challenge to the present-day technology. As an alternative, the 
field can be manipulated by using recently proposed jl^ Nb/AlO^Nb-Al/AlOaj/Nb stacked 
Josephson junctions instead of SQUIDs. In the double Josephson junction device, injection 
current Imj (of the order of several mA) in the control junction induces variation of the 
Josephson critical current Ic (of the order of several /lA). As for a single Josephson junction 
Bj. ~ Ej = /c$°/(27rc) ($° = hc/2e is the flux quantum, h is the Planck's constant, and 
c is the light velocity) manipulation of Ic allows one to control the B^ field. In the experi- 
ment ^3] Ic = (and B^ = 0) for Jj„j ^ 2.5 mA. Using the electrical signals instead of the 
localized magnetic fluxes for conrolling the B^ field would allow one to speed up the quan- 
tum computer operations and to simplify the design. The results presented in this paper are 
applicable also to the quantum computer based on the stacked Josephson junctions (instead 
of SQUIDs) serving as qubits. 
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II. ARRAY OF CAPACITIVELY COUPLED SQUIDS AS A QUANTUM 

REGISTER 

A schematic plot of an array of capacitively coupled SQUIDs is shown in Fig. 1. The 
ith SQUID with i = 0, ■ ■ ■ , N — 1 corresponds to the ith qubit. The SQUIDs have Josephson 
energies Eji and capacitances Cjj. Each SQUID is connected to the control gate voltages 
Vi via a gate capacitor Cgi. The ith intermediate qubit is connected to its two neighboring 
{i ± l)th qubits via the capacitors Cj^j-ti, where Cj^j+i are the off-diagonal elements of the 
capacitance matrix. In Fig. 1 we use the notation Cd = Cj.j+i. The end 0th and (A^ — l)th 
SQUIDs are connected to one, respectively, the 1st and (A^ — 2)th SQUIDs. The Hamiltonian 
corresponding to the charging energy of this system can be written as 

N-l 

Hc(t) = 4 5: [n, - n,,(t)]r,,[n, - n,,(t)], (3) 

i,j=0 

where rij is the operator for the total number of Cooper pairs in the ith SQUID, and ngi{t) = 
CgiVi{t)/2e is the charge (in units of 2e) induced on the i-th qubit by the corresponding dc 
electrode; Tij is related to the inverse of the capacitance matrix Cij of the system by 

r., = ^Cr^- (4) 

To illustrate the structure of the capacitance matrix we assume that all SQUIDs are identical 
with the same capacitances C^j = C^, Cji = Cj, and Cd = Cc- (This assumption is made 
for this example only.) Then, the capacitance matrix is 



Cij — Co 



6ij{l + gr]) - T]J2ki+d 



(5) 



where Cq = C^ + Cj, r] = Cc/Cq; 6ij is the Kronecker delta-symbol; d = ±1 and g = 2 
for the intermediate qubits; d = 1 and g = 1 for the 0th qubit; d = —1 and g = 1 for the 
last (A^ — l)th qubit. If other devices (for measurement) are attached to the SQUIDs, their 
capacitances must be added to Cq. 

Since the matrix Tij is the inverse of a tridiagonal matrix, it has nonzero matrix elements 
on the second, third and other diagonals which characterize the long-range interaction be- 
tween the qubits. However, for <^ 1 the off-diagonal elements of T^j decay exponentially 
as T]!*"-'!, so that the influence of the long-range interaction can be reduced by taking the 
coupling capacitances Cc to be much smaller than the on-site capacitances Cq. (In the 
experiment Q the value of r] was chosen to be r/ ~ 0.05.) 
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Usually, long-range interaction between qubits is difficult to control, and this interaction 
represent a problem for quantum computation. However, in the case when there are only 
two states in the quantum register, we will show that the destructive effect of the long-range 
interaction can be completely suppressed by the proper choice of protocol parameters, so 
that the condition r/ <C 1 is not required in this case. Since the clock speed of the quantum 
computer is proportional to t], increasing the value of rj would allow one to increase the 
speed of implementation of a quantum algorithm. This particular example with two states 
can be important for benchmarking a scalable superconducting quantum computer device. 

In order to obtain the effective Hamiltonian in spin-| notation (see, for example, 
we make the substitution n^j = ?t,° + rrigi and rij = ra^ + nij, where rrij is the operator of 
excess number of Cooper pairs on ith SQUID. The dimensionless charging Hamiltonian (jH)) 
becomes 



C 



Hc(^) = A^cit) = E K - m,,it)n,[m, - m,,(t)], (6) 

i,j=0 

where Uij = rjjCo/(2e^) and Cq is the on-site capacitance for the 0th qubit. It is convenient 
to introduce the spin operator Sf = 1/2 — rrij with projections s| = ±1/2 on the z-axis. We 
treat nij as an operator with eigenvalues and 1. Then, the state with the eigenvalue mj = 
corresponds to the spin state ( ^ , and the state with the eigenvalue rrii = 1 corresponds 



0, 



to the spin state |^ J [l|. Using the relation between rrij (nij) and (Sj), the symmetry 

/I 0\ 

property Uij = Uji, and the fact that (Sf )^ ~ ^ i j ' Hamiltonian of the whole 
system can be written in the form 

N-l N-l N-1 

H(t) = - E Bm^t - E ^rws- + E ^^.s-sf (7) 

i=0 i=0 'J^O 

Here we omitted the term Y.%Zl Uij[mgi{t) - 1/2] [mgj{t) - 1/2] + (1/4) E^Io^ ^u, which does 
not contain the spin operators and does not influence the quantum equations of motion. 
The controllable and fields are expressed through the parameters of the model: 

N-1 p p 

Btit) = E U.A-^ - 2m,,(t)], B^{t) = -^,Ej,m)] = cos [7r$.(t)/$°], (8) 

j=0 ^ 

where E^^ is the Josephson energy of each of the two Josephson junctions of the i-th SQUID, 
^i{t) is the magnetic flux through the i-th SQUID. From the second expression in Eq. (jH} 
it follows that i?f = for non-zero flux $j = {2k - l)$°/2, where k = 1,2,.... 
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A. Two qubits 

Some basic properties of the system can be understood from an exact analysis of the 
SQUID chains containing two or three qubits. For the system of two capacitively coupled 
qubits, the capacitance matrix is 

V -Ce Ci + Cj V ~V a + vJ 
where a = Ci/Cq; Cq = Cgo + Cjq is the on-site capacitance for the 0th qubit; and Ci = 
Cgi + Cji is the on-site capacitance for the 1st qubit. The matrix Uij has the form 

^ a + {l + a)r^\ ^ l + r^j' ^^^^ 
Since the coupling constant Uqi = ///[a + (1 + a)rj\ is positive, the coupling in the system is 
antiferromagnetic. This is also true for larger qubit arrays. 



B. Three qubits 



Assuming that all on-site capacitances and coupling capacitances are the same for all 
SQUIDs (we make this assumption only for this particular example), the capacitance matrix 
for the system of three qubits is 

\ + ^ -ri 



Cij 



/ Co + Cc — Cc 

— Cc Co + 2Cc — C 

-C, Co + c^y 



Cn 



2r] 



V 



The matrix [/,,, has the form 



V 



n2 



1 + 277 + 772 



?7 



\ 

-7] 
l+T]) 

\ 



For the system of three qubits the coupling constants are 



f/ni = U 



1] + 1] 



12 



1] + 7] 

1 + 377 + 77 V 

if 



:ii) 



(12) 



_ nq\ 

l + 4r/ + 3r/2' 1 + 4^ + 3^2- ^ ) 

As mentioned above, the off-diagonal matrix elements Uij decrease approximately as t/'*"-'', 
where r/ <^ 1. The diagonal components Ui^i, which affect the field field Bf in the first 
expression in Eq. (jH)), depend on the position of the qubit in the chain, i. For a system 
consisting of larger number of qubits, the off-diagonal components also depend on i. 
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III. THE Bf OPERATOR FIELD 



It is convenient to express the Hamiltonian (|7j) in terms of the operator field B^(t), 



N-l 



N-1 



i=Q 



i=Q 



N-l 



m) = E {U^: 



l_2m,,(t)-S^ 



} + Uii[l - 2mgi{t)]. 



(14) 



(15) 



Here different terms correspond to different sources of the Bf(t) field in the location of ith 
qubit. The term Uij[l — 2mgj{t)] is the field produced by the application of the voltage to 
the jth qubit, j i; the term Uull — 2mgi(t)] is the field created by applying the voltage 
directly to the ith. qubit; and the term —UijSj, j ^ i, is the field produced by the jth qubit 
in the location of the zth qubit due to the constant interaction between them. 

Since the B|'(t) field is an operator, it generates different actual fields for different 
states. The field B-{p,t) in the location of the ith qubit for the state \p) is defined as 

B^{p,t) = {p\Bt\p) = E {Ui, [l - 2m,,(t) - s^^{p)]} + Uu[l - 2m,,(t)], (16) 



j=0 



where Sj{p) is the eigenvalue of S| for the state \p). For example, for the state |0) = |00 . . . 00) 
one has s|(0) = 1/2 for all j's and 



Af-l 



B^{o,t) = {o\Bm = j:\u. 



j=0 



--2mgj{t) 



Uu[l - 2mgi{t)]. 



(17) 



Below we assume that the voltage Tngj{t) and the control fiux $j(t) [Bj field) are varied 
only on one site with j = I, i.e. mgi{t) = + Mi{t). These parameters are constant for 
all other sites, mgj{t) = m] and ^j(t) = $°(2A; - l)/2. A; = 1, 2, ... for j 7^ /. Then the 
Hamiltonian (|14j) can be written as 



7V-1 



H^t) = H° + 2Mi{t) E UuSt - Bf{t)S: 



i=Q 



where 



Here 



Af-l 

E B■•°s^ 

1=0 



N-l 



B 



E [U., (1 - 2m° - Sj 



i=o 



aAi-2m\ 



0\ 



(18) 



(19) 



(20) 
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FIG. 2: The static field (0) function of tlie qubit number i for two values of r/, nv- = 
(j = 0,l,...,19), iV = 20. 

is the static Bf operator field. The actual static Bf field for a state \p) is 

Bf{p) = mf\p) = E {u^, 

3 = 
37^i 

where the argument p of B^'^{p) indicates the state number. This field is nonuniform even 
for a uniform qubit chain. To illustrate this, we consider the particular example in which 
all qubits are in the same state |0) = |00 . . . 00) and all static voltages have the same values 
m° = 0. Then from Eq. (j21|) one obtains 

1 N~l 

Bf{0) = U,^, + -Y.U.r (22) 

j=0 

In Fig. 2 we plot the static field -Bf'°(0) as a function of the qubit number i for two values 
of 7] and for = 0, j = 0, 1, 19. The number of qubits is = 20. For a two-qubit chain 
with the matrix (jlO|) and a = 1 one can show that i?o'°(0) = i?i'°(0) ^ 1 — 77/2 for 77 ^ 1, 
which is approximately equal to the values of i?Q'°(0) and B^jf_^{{)) in Fig. 2 for rj = 0.1. 
The static field in Fig. 2 is nonuniform near the edges of the qubit chain despite the fact 
that all applied voltages have the same values = and all qubits are in the same state 
|0,). 




r|=0.5 



6^ ^ 

Ch -o — o- -o — o— o- -o — o- -o —0—0- -0—0- -o — o— cr 



2m^ 



(21) 
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In dimensional units the total static voltage applied to ith qubit is = (2e/Cc,)(n°+m^). 
If = m° = 0, then f° = 0. Due to Eq. © the relation between the dimensional magnetic 
field i?j^'°(0) and the effective static voltage, corresponding to this field, can be written as 
(2eVCo)5,^'°(0) = evf^. For example, the value Bf{0) ^ 0.92 for i = l,2,...,iV - 2 
and 7] = 0.1 in Fig. 2 for Co ~ 500 aF Q] corresponds to the effective voltage v"^^^ = 
0.92 X 2e/Co ^ 0.6 mV. 



IV. THE INTERACTION REPRESENTATION 



We decompose the wave function into the basis states \p) of the unperturbed Hamiltonian 
HO; 

m = E cp(t)b) = E A,{t)e-''^^'\p), c,{t) = A,{t)e-''^'^\ (23) 

r=0 p=0 

where the Planck constant is h = 1, \p) = |nAr_ir?,jv-2 ■ ■ - Ui . . . riiUQ), rii = 0, 1, and 

N-l 

E, = (p|H» = -Y: Bf{p)st{p). (24) 
The dimensionless time t in Eq. is expressed in terms of the dimensional time t as 

For Co = 500 aF, the dimensionless time t = 1 corresponds to the dimensional time t ~ 1 
ps, and one dimensionless energy unit corresponds to 0.64 meV. 

In the stationary field Bl'^{p) and when Bf = 0, the coefficients Apit) do not change 
when time t changes from t' to t", i.e. Ap{t") = Ap{t'), while the coefficients Cp{t) evolve as 
Cp{t") = Cp{t')e~^^p^^"~^'\ The representation of the wave function in which the coefficients 
Ap{t) are used is called the "interaction representation". Since the coefficients Ap(t) are not 
changed during the free evolution of the system, the interaction representation allows one to 
exclude from consideration the dynamics associated with the evolution of the phase of the 
wave function when no pulses are applied. 



V. ONE-QUBIT FLIP 

Let us discuss how to implement a resonant one-qubit rotation of the Ith qubit for the 
state \p) in our computer. Initially at the time to the 2;-component of the magnetic field is 
Bi'^{p), and the Bf field is switched off, B'f{tQ) = 0. The flip is implemented in three steps. 
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(a) One changes the effective voltage Mi{t) (below called voltage) applied to the Ith qubit, 
so that at the end of the pulse of duration ti = ti — the magnitude of M;(t) becomes 
Mi{ti) = Ail. If the transition is resonant Aii is defined by the resonant condition. [See 
Eq. ()Hfi|l below.] From Eq. ()18|1 one can see that after this pulse the field in the location 
of the Ith qubit is 

Blip,t,)=Bf'\p)-2MiUu. (26) 

Since during the time of application of the first pulse the Bf field is turned off, the form 
of the first pulse is not important, and the relevant parameter is the area fl^ Mi{t)dt of 
this pulse which determines phase acquired by the wave function during this pulse, (b) 
One flips the Ith qubit by a rectangular vr-pulse with the amplitude Bl and time duration 
T2 = t2-ti = n/\Bl\. The B[ field is not changed Bf{p,t2) = Bf{p,ti) and Mi{t2) = 
Mi{ti) = Ail. During this pulse the spin flips through the angle vr. (c) One restores the 
voltage applied to the Ith qubit from Mi{t2) = M.i to its original value M/(t3) = during 
the time = — t2. 

First we will discuss the dynamics of the quantum computer during implementation of 
the steps (a)-(c) in terms of the coefficients Cp(t), then we will formulate the result in the 
interaction representation in terms of the coefficients Ap{t). Let the Ith spin of a state \p) 
at the initial time to be in the state and the state |g) be related to the state \p) by a flip 
of the Ith spin. 

After the first pulse one has 

Cp(ti) = e-'^''"^-'^'(P)-'^'cp(to), c,{h) = e-'^''^'-''l^P^+'^"cg{to). (27) 

The total phase (for the states \p) and |g)) Oji^p) and the phase ipj are defined by the 
time-dependent component of the Bf field in Eq. ()18|1 as 



el{p) = 2 / Mi{t)dt Y^Uustip), vl= / Mi{t)dt 



1=0 



Uij. 



(28) 



At the end of the first pulse Mi{ti) = Mi. 

After the third pulse, the probability amplitudes become 

Cp{ts) = e-'''''^'-''f^P^-'^'Cpit2), c,{h) = e-'^^^^-^^?(^')+''^'c,(t2). 
The phases are 



Mi{t)dt 



N-l 



E UhsI{p), 



i=0 



t2 



Mi{t)dt 



Uu. 



(29) 



(30) 
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At the beginning of this pulse, one has Mi{t2) = M.u and at the end of the pulse Mi{t^) = 0. 

The action of the second rectangular pulse is described by the following Schrodinger 
equation: 



where 



icp{t) = [Ep + blip) + MiUa] Cp{t) - -j-c,{t) 
icq{t) = [Eg + blip) - MiUii] c,{t) - -j-cp{t), 



N-1 

bi{p)=2MiY.Uusl{p). 

i=0 



(31) 



(32) 



In order to make the transition \p) — > \q) resonant, the diagonal elements in Eq. ()31|) 
must be equal to each other. The condition 



Eg-Ep = 2M\Uii 



(33) 



defines the resonant value of the applied voltage M.\. We now calculate the energy difference 
Eq — Ep using Eq. Since the Bf'^{p) = Bf'^{q) and sf{p) = —sf{q), the contribution 

to the energy difference directly related to the flip of the Ith qubit is equal to Bf'^{p). The 
states of other qubits with i ^ I are the same for both states, sf(p) = sf(g), but the fields 
given by Eq. (j^D), are different due to the influence of the Zth qubit measured by the 
matrix elements Vu^ 



N-1 



j=0 



Finally, we obtain 



N-1 



Eq - Ep = Bf{j>) -Y.Uus\{p). 
From Eqs. and (jH^j) the resonant value of the applied voltage is 



(34) 



(35) 



2Un 



Bf{p) - E UhsUp) 



N-1 



i=0 



(36) 



Note, that the resonant value of the B[ field is not zero because of the constant inter-qubit 
interaction. From Eqs. (j26j) and (j36|) it is equal to 

N-1 



Bt{pM) = Bf{p)-2MlUii 



1 



mi u 



E Uustip). 



(37) 
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To find the solution generated by the second pulse, introduce the new coefficients Dpij) 
and Dq^r) as 

Cp{t) = Z}p(r)e-'[^''+'"(*')+^'^"l", 

Cg{t) = Z}^(r)e-'[^^+^'(f)--^'^"]", (38) 
where t = t — ti. Then, for the coefficients Dpi^r) and Dg^r) one obtains equations 



lDp{T) 



Dp{r)e^^'\ 



where 



N-l 



Ai = Eq-Ep- 2MiUu = Bt'\p) - J2 UusUp) - 2MiUu. 



(39) 



(40) 



i=0 



The condition A/ = is satisfied if the pulse is resonant [see Eq. (jH^ ]. 

If initially the Zth qubit is in the state 0, -Dp(O) = Cp(ti) and -0^(0) = 0, the solution of 
Eq. dSni) is 

Dp{T2) = Cpiti) 



(XlT2\ , .A^ . (X1T2 



Dq{T2) = Cp (ti)'i--^ sin 



^Bl ^.J^ljA 



(41) 



where Xi = \J (Bfy + (A;)^. From Eq. (j4T| one can see that if A; = and T2 = 7r/\Bf\ there 
is a complete transition between the states \p) and Below we call A/ for state \p) the 
detuning for this state. 

We now express these results in the interaction representation in terms of the coefficients 
Ap{t). Let the state \p) with Zth qubit in the state |0/) be populated at time t = to. From 
the second expression in Eq. (|^ and the first expression in Eq. ^I7\ one obtains 



Cp(to) = Ap{to)e 



-iEpto 



(ti) = Ap{to)e-'^'^''-'^'^P^-'^' . 



(42) 



Taking into consideration Eqs. (|HH|) . (HOJ, and (jH}, one has 



Cp{t2) = Ap(to) 



( XiT2\ , .A, . / X1T2 



■{ 



i\ -Ept2-eiip)-y,i- 



^ A(+\-B} . (hr2\ i{-E,t2-e]{p)-^]- 
Cg{t2) = Ap{to)i—sm I I e 1^ 



bi{p)+MiUa + ^ 



T2 + {^l+2MlUu)t2^ 



(43) 
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Finally, using Eq. ()29|) one obtains in the interaction representation 



( \lT2\ , Al . ( \lT2 



,i[eKp)-ef{p)-¥'f] 



A,(t3) = Ap(to)^-r^sin 

A; \ z, / 

where the common phase 
depends on the parameters of the first two pulses and on the initial state 



bi{p)+MiUu + Y 



T2 



(44) 



(45) 



A. Calculation of the detuning A/ 

Assume that initially there are two states, |r) and in the quantum register and 
assume = for the state \R) and A; 7^ for the state |r): the pulse is resonant for 
the state \R) and nonresonant for the state |r). The value of A; for the state |r) can be 
found from Eqs. (^^1 and (jlO)), 



N-l 



Ai{r,R) = -2Y.Uu[st{T) - Sim- 



(46) 



i=0 



This equation can serve as a general definition of the detuning A^ for an arbitrary state |r) 
expressed through the eigenvalues of S| for this state and the eigenvalues of this operators 
for the state for which the transition is resonant. From Eq. ()4fij) one can see that 
in general A/(r, i?) is independent of the static voltages m^, and for a resonant transition 
A;(i?,i?) = 0. 



VI. PROTOCOL FOR CREATION OF ENTANGLED STATE 

Let the initial state be the state |0Ar_i0Ar_2 • • • OiOq). If ^ 1 this is the ground state 
when the applied field is oriented in the positive 2;-direction, i.e. when < m° < 1/2. 
(See Fig. 2 for the case m° = 0.) Using the Hadamard transform Gh we split the ground 
state into two states 

GH|0iV-l0Ar-2 . . . OiOo) = ^{\^N-l^N-2 ■ ■ ■ OiOq) + |0;v-l07V-2 • • • OiIq)). (47) 
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Here and below we omit the total phase factor. During the next step, we flip the first qubit 
in the excited state and do not flip the same qubit in the ground state to obtain the state 

-^(|0jv_i0jv-2 • • • OiOo) + |0;v-i07v-2 • • • Oslilo))- (48) 

Repeating the latter procedure for the remaining N — 2 qubits we will obtain the entangled 
state: 

-^(|0jv^i0jv-2 • • • OlOo) + \In-i1n-2 . . ■ lilo))- (49) 
We now define the parameters of the pulses required to implement this protocol. 



A. Hadamard gate 

Assume that initially, at time to; there is only the ground state |0) = |00. . .00) in the 
register. The Hadamard transform is implemented using the sequence of the pulses (a)-(c) 
described in SecjVl In step (b) instead of vr pulse we apply a 7r/2-pulse. 

The form and duration of the first pulse are not important, since they affect only the 
total phase of the wave function. What is important is the value of the i?o(0,ti) field after 
this pulse. To make this field equal to the resonant value (resonant pulse) one should apply 
the voltage 



Ml 



N-l 



4 = 1 



(50) 



2Uoo 

where s^(0) = 1/2, and we used Eq. (jSHl)- During the second pulse one keeps the gate voltage 
constant. For the case m° = one has Ai^ = 1/2. The second rectangular 7r/2-pulse with 
the amplitude Bf has the duration r2 = 7r/(2|i?/|). A specific value of Bl is not important. 
The only condition is Bf ^ 0. 

During the third pulse the voltage is changed from A4q in Eq. (j50|) to its original value 
which is equal to zero. To equalize the phases of two states of the superposition, the condition 

-^1 = ^1 + [Ao(0, 0) + 2MlUoo]t2 + ^ (51) 

must be satisfied [see Eq. (jSj)]. Here Ao(0, 0) = 0. If we set initially to = 0, then 

t2 = ri + r2 = ri + ^. (52) 

This gives us the last parameter of the Hadamard gate 

= -MoUoo U + ^) - ^. (53) 
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B. Conditional flip 



All other steps of the protocol implement flips of the 1st, 2nd . . ., (A^ — 2)th, (A^ — l)th 
qubits in the excited state and suppress flips of these qubits in the ground state. We 
will derive the parameters required to implement this operation for the flrst qubit. The 
parameters required to flip the other N — 2 qubits can be obtained in a similar way. Note 
that the parameters required to flip different qubits are different even for a homogeneous 
spin chain because the static fleld B^'^{R) is different for different qubits i and states 
\R), and because the phase factor exp[i2J^iUiit2] in Eq. (jUj) depends on the history of the 
excited state. 

Before the pulses (a)-(c) are applied, there are two states in the register |r) = |00 . . . 00) 
and l-R) = |00 . . . 01). As for the Hadamard transform, the form and duration of the first 
pulse is unimportant. The value of the voltage after this pulse is 



2Uu 



N-1 

Bf{R) - y: uusm 



i=0 



(54) 



If = 0, then 

The second rectangular vr-pulse with the amplitude B} has the duration T2 = n/{\Bl\). 



(55) 



In order to suppress the nonresonant transition |00. 
state, the value of \Bl\ must satisfy the 27rA'-condition 

Ai (r,i?) 



..000) 
3n I20I 



K=h2, 



|00 . . . 010) from the ground 



(56) 



~V4:K'^ - 1' 

where Bl can be positive or negative. For this particular pulse, using Eq. ()46|1 . one obtains 
Ai (r,i?) = -2f/io. If Bl satisfies Eq. §^ the value of the sine in Eq. ()44|) becomes zero, 
and the nonresonant transition is suppressed. 

The voltage Mi{t) is switched off, A4{ 0, and the phase is corrected by the third pulse, 
for which the condition 



ei(r) 



N-1 



<^ll i=0 



-Ktt 



2 

1 - TT E Uusm 



IT 



+ 0i(i?)+2A4if/nt2+2 
(57) 
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must be satisfied. Here t2 — Ti + ti + r2, and Ti is the time of the beginning of the current 
three-pulse sequence. Finally, the phase is 

3 _ ei(r) - Qi{R) - 2MiUut2 + {K- l)f _ Qi{r) - QijK) - 2MiUut2 + {K ~ l)f 

(58) 



VII. COMMENTS ON PRACTICAL IMPLEMENTATION 

The protocol discussed in this paper provides the parameters for exact implementation of 
entanglement in a chain of many coupled SQUIDs. We demonstrated that the phase of the 
wave function can be controlled by controlling the gate voltages. In a practical situation one 
of the objectives can be the creation of the entangled state with a relatively large constant of 
interaction between qubits (but not using this state for implementation of more complicated 
quantum algorithms). If the creation of the entangled state is the purpose of an experiment, 
then the phases of the entangled states can not be important. In this situation, there is 
no need to control the phase (pf, and the important parameters are: (i) the value of the 
gate voltage A4i apphed to the qubit to be flipped, (ii) the amphtude Bl and (iii) the 
time-duration T2 of pulse [controlled by the magnetic flux ^i(t)]. 

We now make some remarks concerning the specific field configuration that we chose 
in this paper for implementation of quantum logic. First, all qubits are initially placed in 
a permanent static B^ field, Bj'^ ^ created by static voltages m° and the inter-qubit 
interaction. Before a flip of the Ith qubit by a B^ pulse, the voltage apphed to this qubit is 
switched to the resonant value, and after the S^-pulse this voltage is switched to its initial 
value. As follows from the results of this paper, the nonzero static field introduces an 
additional parameter, or degree of freedom: by changing this field it is possible to control 
the phase of the wave function. 

The qubit rotations are implemented by the rectangular -pulses because during imple- 
mentation of these pulses the effective Bf field in the location of the Ith qubit is not equal 
to zero for some quantum states of the superposition. That is why the result of the action 
of a S^-pulse on the states with nonzero Bf field depends on the form of this pulse. The 
analytic solution is known for a rectangular pulse, so the rectangular form of the i?f -pulse 
is the simplest possible choice that allows one to analyze the basic properties of the system 
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and optimize the parameters of the pulses. 

The interaction representation used in this paper for modehng of quantum logic is em- 
ployed because it effectively eliminates from consideration free quantum dynamics of the 
quantum states of the superposition. Namely, since different states of the superposition 
have different energies Ep, they acquire different phase factors even in the case when the 
applied fields are not changed. The interaction representation allows one to choose the static 
voltages m° arbitrarily. In particular, one can choose m° = 0. (See, for example. Fig. 2 
of this paper.) If n° = 0, this would allow one to decrease the decoherence effect of the 
gate because the gate voltage is switched on and off only during the logic operation on a 
specific qubit. Also the time between the pulses can be chosen arbitrarily because the wave 
amplitudes in the interaction representation are not changed during this time. 

VIII. LONG-RANGE INTERACTION 

For two states in the register, the system with long-range interaction between qubits 
is exactly solvable. However, the universal quantum computation involving an arbitrary 
possible number of states in the quantum register is supposed to be implemented by taking 
into consideration only nearest neighbor interaction. The long-range interaction in the 
general case results in generation of error states. These errors are small only when the value 
of 7], characterizing the coupling strength, is small. For example, the value t] = 0.1 can 
be considered as large because the long-range interaction produces the error of the order 
?7 = 0.1 in the probability amplitude and phase of the wave function. In order to show this, 
consider Eq. (jUJ. Due to Eq. (jlUj) the contribution of (/ — 2)th and (/ + 2)th qubits to the 
detuning A is of order r/^; the contribution of (/ — 3)th and (/ + 3)th qubits to the detuning 
A is of order rj^, and so on. If the parameters of the pulses are calculated by taking into 
consideration only nearest neighbors, instead of the resonant transition with A = one has 
an almost resonant transition (for 7] <^ 1) with the nonzero detuning A ~ created by 
(/ — 2)th and (/ + 2)th qubits. Since in Eq. (jUj) A; ~ i], the relation Ai/\i is of order of 77; 
the value of the sine (for almost resonant transition) is of order of unity; so that the error 
in the probability amplitude and phase of the wave function is of order rj. 

In order to calculate the error introduced by the long-range interaction, we numerically 
simulated the creation of entanglement. The parameters of the pulses were calculated by 
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Pulse number 



FIG. 3: The probability error as a function of pulse number. N = 100, K = 1; m?- = 0, 
i = 0,l,...,A^-l. 

taking into consideration only F nearest neighbors, F = 1, 2, . . .. We define the probabihty 
error as 



\Mt)?-\ 



+ 



\Mt)?-\ 



(59) 



and plot it as a function of the number of pulses in Fig. 3. The total number of 
pulses is equal to the number of qubits in the chain. In Eq. ()59|) |r) is the ground state 
and \R) is the excited state. 

The probability error V{t) in Fig. 3 generated by one pulse is equal, on average, to 
2.87 X 10-^ (9.77 x 10"^) for F = 5 and r/ = 0.5; and V = 2.93 x lO^^ (1.96 x lO^^) for 
F = 1 and rj = 0.0014. In brackets we indicate the values calculated using the analytical 
estimates presented in the beginning of this section. Namely, the error in the probability 
amplitude is of the order rj^, and the error in the probability is of order r/^^. The poor 
correspondence between estimated and numerical results for t] = 0.5 follows from the fact 
that for relatively large rj the off-diagonal matrix elements in the matrix Uij decrease faster 
than T]!*"-'!, i ^ j. 

As follows from Fig. 3, in order to decrease the error one can (i) take into consideration 
a larger number of nearest neighbors or (ii) decrease the coupling constant rj. The results 
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0.008 



0.006 

^ max 
0.004 

0.002 



20 40 60 80 100 

Number of pulses 

FIG. 4: The maximum value of rj (y-axis) required to implement the entanglement protocol using 
N pulses (j;-axis), with the probability error less than Vq = 10"^. F = 1, = 1; m° = 0, 
j = 0,1, . . . , N — 1. The solid curve is obtained using Eq. H44|) and the dashed curve is obtained 
using Eq.(|HU. 

presented in Fig. 3 can be used to estimate the probability error in more complex protocols 
with the same values of F and t], since one can assume that this error mostly depends on 
the number of pulses and not on a particular protocol. 

As shown in this paper, the effects of the long-range interaction can be completely com- 
pensated by choosing the optimal values of the and B^ fields when the number of useful 
states in the quantum register is equal to 2. In the general case of an arbitrary (possible) 
number of states in the register, only the nearest neighbors can be taken into consideration, 
and the long-range interaction can be suppressed only by decreasing the value of rj. In 
Fig. 4 we plot the maximum value ?7max required to implement the protocol for creation of 
entanglement using N pulses with the probability error less than Vq = 10~^ when only the 
nearest neighbors are taken into consideration, F = 1. 

The value of ?7max in Fig- 4 can be estimated analytically. As follows from Fig. 3 the error 
accumulates linearly with the number of pulses, so that we can write 

Vo = NV = ATr^Lx, (60) 



numerical 

analytical 
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where V is the error generated by one pulse. From Eq. ()6U|) the maximum value of rj required 
to implement an A^-pulse protocol with the accuracy Vq = 10""^ is 



From Fig. 4 one can see that our analytical estimate for r/max is close to the exact result 
obtained using numerical simulations. 

Decreasing rj decreases the error. On the other hand, since rj defines the clock speed 
of the quantum computer (the time of implementation of one pulse is of order of vr/r^), 
decreasing rj slows the computer down. This can lead to an accumulation of errors introduced 
by the environment. The optimal value of f] can be estimated by taking into consideration 
both the effect of the long-range interaction and the influence of the environment. 



We considered in this paper an optimal implementation of quantum logic operations 
for a scalable superconducting quantum computer with constant inter-qubit interaction. 
The protocol for creation of entanglement with arbitrary number of qubits was analyzed in 
detail, for a relatively large interaction constant. A possible application of our approach 
for implementing universal quantum logic for more complex algorithms by decreasing the 
coupling constant and, correspondingly, decreasing the clock speed was discussed. The 
errors introduced by the long-range interaction for the universal logic gates are estimated 
analytically and calculated numerically. The coherent charge oscillations in this model have 
already been observed experimentally |^J. A further feasible accomplishment would be the 
experimental creation of entanglement in a system with three qubits. A demonstration of the 
entanglement in the potentially scalable QC architecture with the superconducting qubits 
would be an important step toward the experimental implementation of a scalable QC with 
many qubits. 




(61) 
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APPENDIX A: CREATION OF ENTANGLEMENT WITH TWO QUBITS 

Let us calculate the parameters required to create the entanglement in the chain of 
two qubits starting from the state |00) and using the parameters from the experimental 
system j^. We will not calculate the phase iff required to implement the phase correction, 
since it depends on particular forms of the gate pulses Mi{t), and calculate the parameters 
required to create the entanglement without phase correction. There are three devices 
attached to ith {i = 0, 1) SQUID: superconducting electrode with the capacitance Cji, 
electrostatic gate with the capacitance Cgi, and a probe with the capacitance Cbi- We 
assume that during implementation of the protocol the probe voltages Vf,i are switched off, 
so that the charges Qbi = CbiVbi induced by the probes are equal to zero. The capacitances 
are 

Cjo = 620 aF, Cgo = 0.60 aF, 0^ = 41 aF, Cq = Cjo + Cgo + Cbo = 661.6 aF. (Al) 

Cji = 460 aF, Cgi = 0.61 aF, Cbi = 50 aF, Ci = 510.61 aF. (A2) 
The coupling capacitance Cc and the coupling constant rj are 

C 34 

a = 34 aF, r/ = = — — ^ 5.139 x 10"^ (A3) 

Lo 661.6 

The coupling matrix Uij, i,j = 0,l, is given by Eq. (fTUj) . where 

C, _ 5ia« ^ 0.7718. (A4) 

Co 661.6 ^ ^ 

The dimensionless parameters rj and a and the dimensionless static voltages tUq and 
completely define the parameters of the protocol. 

The value of the voltage Aio applied to the zeroth qubit is given by Eq. (jHEj), 

A<5 = i - m° - ^m?. (A5) 



The amplitude Bq of the applied field is arbitrary, and the time-duration of this pulse is 

TT + 2'n-k 



T2= ^,^,1 , A; = 0,l,2,.... (A6) 
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The gate voltage Aii for the second pulse is 

M[ = l-m'i + ^{l-0- (A7) 

The magnitude of Ai calculated using Eq. is —2Uoi, and the amplitude of field is 

Bl= , K=l,2,.... (A8) 



The time-duration of the second B^ pulse is 

TT + 27r/c 



r2= 1^11 , A; = 0,1,2,.... (A9) 



APPENDIX B: CREATION OF ENTANGLEMENT WITH THREE QUBITS 

The matrix Uij for the system of three identical qubits is given by Eq. ()12|) . The protocol 
presented in this Appendix is also valid for a system with nonidentical qubits. The gate 
voltage Aio for the first B^ pulse is 

The magnitude of B^ is arbitrary, and T2 is given by Eq. ()A6|) . 
The gate voltage TVli for the second pulse is 

Mi = i-m? + ^(l-m°)-^m°. (B2) 
/ Uii Uu 

The magnitudes of Ai = —2Uoi, Bl, and T2 for the second pulse are given by the same 
formulas ()A8|) and ()A9|) as for the case of two qubits. 
The gate voltage for the third pulse is 

Ml = l-mt + ^{l-m',) + ^{l- m% (B3) 

^ U22 '-^22 

The detuning A2 for the third pulse is —2U20 — 2U21, and the amplitude of B^ field is 

'~ V^K^'' ^-1'2,..., (B4) 



The time-duration of the third B^ pulse is 

vr + 27r/c 
LB 



T-2= ,^11 , fc = 0,l,2,.... (B5) 
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